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Intrinsic instability of sonic white holes 
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Artificial black holes, such as sonic holes in Bose-Einstein condensates, may give insights into the 
role of the physics at the event horizon beyond the Planck scale. We show that sonic white holes 
give rise to a discrete spectrum of instabilities that is insensitive to the analogue of trans-Planckian 
physics for Bose-Einstein condensates. 
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Hawking predicted jl) that black holes generate ther¬ 
mal radiation due to quantum effects [§| [|. Hawking’s 
effect bridges two vastly different areas of the physi¬ 
cal sciences — cosmology and quantum mechanics, but 
unfortunately, the effect is too feeble to be observable 
for the known solar-mass or larger black holes. How¬ 
ever, condensed-matter or optical analogs may be able 
to demonstrate the equivalent of Hawking radiation in 
the laboratory @, g, g, | g, g, [0, 0 0. 0.00.0. 
In understanding the mechanism of laboratory analogs 
one could perhaps gain insight into the anatomy of gen¬ 
uine quantum black holes. For example, Hawking’s ef¬ 
fect appears to rely on a theory that predicts its own 
demise [0] - radiation outgoing from the event horizon 

seems to originate from wavelengths beyond the Planck 
scale where the physics is unknown. Yet the theory of 
condensed-matter analogues indicates that Hawking’s ef¬ 
fect is robust against trans-Planckian physics J0] . In this 
Letter we show that white holes with effective surface 
gravity a give rise to a discrete spectrum of dynamical 
instabilities with decay constants 

7 = 2 na , n € Z , (1) 

whereas black holes may be stable. According to Cor¬ 
ley and Jacobson |J, a black-hole white-hole pair should 
act as a black-hole laser, i.e. as an amplifying medium 
for Hawking radiation with the two horizons forming a 
resonator, giving rise to a dynamic instability of Hawk¬ 
ing radiation. This instability has been attributed to 
physics beyond the Planck scale. Here we point out that 
the discrete spectrum of a single white hole is insensi¬ 
tive to the equivalent of trans-Planckian physics within 
the model we have employed, the sonic hole in a Bose- 
Einstein condensate. Such a model may be regarded as 
the drosophila of the artificial black holes 0, [0], be¬ 
cause it is the simplest system to study theoretically and 
it is within experimental reach. 

Sonic holes are inspiring, because of the mathematical 
equivalence between the propagation of sound in fluids 


and of scalar waves in general relativity g, |. Consider 
sound waves in an irrotational fluid of density po and flow 
u. The velocity potential ip and the density perturbations 
p s of sound obey the linearized equation of continuity and 
the linearized Bernoulli equation |l9| 

dtps + V • (up s + poSPp) = 0, (2) 

(d t + u • V)v? + c 2 — = 0. (3) 

Po 

The resulting wave equation in relativistic notation reads 

D v D v ip = - ? L= = 0 (4) 

with the effective space-time metric [g, in d spatial di¬ 
mensions, 



Nonuniform profiles of po, c 2 and u may generate effective 
space-time geometries that are sufficiently rich to possess 
event horizons g. At a sonic horizon the flow exceeds 
the speed of sound. Sound waves propagating against 
the current freeze and, in turn, their wavelenghts shrink 
dramatically. We are interested in effects dominated by 
this wave catastrophe J0. In this case we can use the 
simple one-dimensional model 

u = —c + az . (6) 

Here £ denotes the spatial coordinate orthogonal to the 
horizon at z = 0, a characterizes the surface gravity or, 
in our acoustic analog, the gradient of the transonic flow, 
and po and c are assumed to be constant. Strictly speak¬ 
ing, we should complement the flow component (g) in 
the z direction by appropriate components in the x and 
y directions, in order to obey the continuity of the flow. 
But as long as we focus on effects on length scales smaller 



2 


than \c/a\ we can ignore the other dimensions of the fluid. 
Equations (j|) and (|ij) have the solutions 


( P = Vq{t), p s 


Po dip 
caz dr ’ 


T= M (7) 

a 


with the arbitrary function tp^ and the constant z^, 
describing wavepackets propagating against the current. 
Such wavepackets are confined to either z < 0 or z > 0, 
depending on the sign of z^, which indicates that the 
place where the flow exceeds the speed of sound, z = 0, 
indeed establishes the acoustic equivalent of the event 
horizon. Depending on the sign of a, two cases emerge: 
the sonic black hole and the white hole |nj. The black 
hole is characterized by a positive velocity gradient a , 
the flow goes from subsonic to supersonic velocity at the 
horizon, whereas a is negative for the white hole |l(J 
where the flow slows down from supersonic to subsonic 
speed. No sound wave can leave the supersonic zone of 
the black hole and no sound can enter the white hole. 

Experimental tests of the subtle quantum effects of 
the acoustic horizon ||, will take the best superflu¬ 
ids currently available, Bose-Einstein condensates of di- 
pi^ . Moreover, a condensate is one of the 
simplest physical systems to draw theoretical analogies 
between the quantum mechanics of fluids and quantum 
effects of gravit y |j| . The mean-field wavefunction tpo of 
the condensate ] 22 | represents an irrotational fluid, 

V’o = y/Po e * So > u=-VS 0 , ( 8 ) 

m 
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where m denotes the atomic mass. Sound waves are per¬ 
turbations ip — V’o °f the fluid, with 


Tp = \/ r pe lS , p = po + p s , S = S 0 + —ip. (9) 


Sound quanta, phonons, are elementary excitations char¬ 
acterized by the Bogoliubov modes u n and v n [p2| p3| , 
with 


The Bogoliubov modes differ on the two sides of the hori¬ 
zon, depending on the analyticity of z lu ^ a , which reflects 
the independence of the two sides of the sonic horizon. 

Transonic flows are notoriously plagued by hydrody¬ 
namic instabilities ||l9[ , unless they are generated in ap¬ 
propriately designed nozzles such as the Laval nozzle |2{| 
of a rocket engine. Consider unstable elementary excita¬ 
tions corresponding to Bogoliubov modes with complex 
frequencies u> j24| , quasinormal modes p 6 | . The dynamic 
equations of the excitations, the Bogoliubov-deGennes 
equations |2^, p3| |24|, possess a four-fold symmetry in 
the complex frequency plane [p"ol 241: If to is the fre¬ 


quency of a solution then there exist solutions for the 
frequencies and —uj* as well. The Hermiticity of 

the underlying many-body Hamiltonian causes this sym¬ 
metry ju, 241. For unstable modes the Bogoliubov scalar 
products | 2 :i|. | 2 ~| must vanish, 


/ -t-oo 

{u* n u n ' - v*v n >) dz = 0, (13) 

-oo 

/ +oo 

{v n U n > - UnVn') dz = 0, (14) 

-oo 


because, in the Bogoliubov-deGennes dynamics [£2j |2fj, 
p~| the scalar products are stationary |23|, p4|, whereas 
the modes © are growing or decaying for complex fre¬ 
quencies. Unstable modes in other field theories are sub¬ 
ject to similar requirements [p 6 |. As a consequence of the 
analytic property ( 12 ), condition (jl^) is satisfied in the 
case of purely imaginary frequencies, 


ui = * 7 . (15) 

Modes ( |lT| ) with positive 7 /a are localized near the hori¬ 
zon at z = 0. We deform the contour of integral © to 
a large semicircle with radius r around the origin on, say, 
the lower half plane, and get 


ip = ip 0 + e lSo (u n + v* n ). ( 10 ) 

We determine the modes from the solution (|) of the 
hydrodynamic equations (|],[|) by comparing the repre¬ 
sentations (||) and © in the limit of \p 8 /po\ 1 , 

IvV'S'ol h/m ©. We find the single-frequency Bo¬ 
goliubov modes of the transonic flow 

= ^ (11 > 

Consider the analytic continuation of z lul ^ a to complex z. 
Suppose that z lul / a is analytic on either the upper (+) or 
lower (—) half plane. Consequently, we get for real and 
positive z 

= g— K(±u>/ot) z iu)/a ' n 2) 



~ W£(7^)/ V)^/“*. ( 16 ) 


The integral vanishes for the decay constants (jlj). Be¬ 
cause of the frequency symmetry of the Bogoliubov- 
deGennes equations 0-0 solutions for negative 7 /a 
must exist as well, although we cannot represent them 
as the acoustic Bogoliubov modes ©. because they 
would grow in space. Consequently, intrinsic instabili¬ 
ties of sonic horizons, if any, correspond to the discrete 
spectrum (Q). Let us scrutinize the assumptions made. 

Close to the horizon the wavelength of sound would 
shrink beyond all scales and the density ( 0 ) would di¬ 
verge, if the wave equation (Q) were universally valid. In 
the short-wavelength limit we can describe the Bogoli- 
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ubov modes in the WKB approximation 00- 


ti„ = U n exp (i / kdz — iut I , 


v n = Vn exp * / kdz — iut 


(17) 


(18) 


We obtain the wavenumber k from Bogoliubov’s disper¬ 
sion relation |2^| in moving condensates, taking into ac¬ 
count the Doppler effect, 

(u - uk) 2 = c 2 k 2 ^1 + , k c = ^ . (19) 


The group velocity 


du 

dk 


u + v , 



( 20 ) 


indicates that the acoustic Compton wavenumber k c de¬ 
fines the trans-acoustic scale beyond which the excitation 
velocity in the fluid, v', deviates significantly from the 
speed of sound. Consider the turning point Zq where the 
group velocity (|2^) vanishes. For single-frequency modes, 
the excitation flux is conserved |24|, [27| , 

d z (U 2 -V 2 )v = 0. (21) 


Consequently, the amplitudes U n and V n diverge at the 
turning point. If the trans-acoustic scale k c were zero the 
horizon itself would be the turning point. Therefore, we 
use |zo I to estimate the spatial range of the trans-acoustic 
region. For elementary excitations [^2| pd[, 


hu 
me 2 


( 22 ) 


is a small parameter. We expand Zo in a power series in 
e 1 / 3 and solve v = 0 to leading order Q, 


zo 


lityzi(i) 

a2 V 2/ 


\ 2/3 


(23) 


The trans-acoustic region \z\ < \zq\ is small compared 
with |c/a|. Therefore, neglecting this region in the inte¬ 
gral (|fj|) is justified. Consider the solutions of the disper¬ 
sion relation ( |l9| ) in the acoustic regime and in the trans¬ 
acoustic extreme. We obtain four branches of this fourth- 
order equation that we characterize by their asymptotics 
on one side of the horizon. As long as k 2 is much smaller 
than k 2 we get the acoustic asymptotics 


k 


u 

k±c ’ 


and in particular, 



az 


(24) 


(25) 



FIG. 1: Stokes lines (dotted lines) of unstable elementary ex¬ 
citations at a sonic white hole. The excitations have purely 
imaginary frequencies. The Stokes lines originate from the 
turning points (dots). The jagged lines indicate the branch 
cuts of the inverted dispersion relation k = k{u). Three 
branches are connected by the branch cuts and may be con¬ 
verted into each other at the Stokes lines, the acoustic branch 

(A) with asymptotics (25) and two trans-acoustic branches 

(B) and (C) with asymptotics (26). The pairs of letters indi¬ 
cate which branches of the superposition are potentially con¬ 
verted into each other. The first letter of each pair identifies 
the exponentially dominant branch (determined numerically). 
The picture shows that we can construct Bogoliubov modes 
such that the (B) and (C) branches are not present on the 
lower half plane, without causing contradictions. Therefore, 
the unstable elementary excitations of the sonic white hole do 
indeed obey the asymptotics (25) on the lower half plane. 


for sound waves propagating against the current. In the 
other extreme, when k 2 is much larger than k 2 , one finds 

0 

k ~ ±2k c y/u 2 /e 2 - 1 + UU 9 . (26) 

c z — u- 

Our analysis of the instabilities is justified if the Bogoli¬ 
ubov modes obey the asymptotics (|25|) in one of the com¬ 
plex half planes. 

At a turning point zo two wavenumber branches of the 
dispersion relation (0) coincide |^, |24j. Therefore, a 
solution that starts from a particular k\ branch on one 
side of zo may become converted into a superposition 
of the two modes k\ and &2 with the common turning 
point. In other words, turning points may cause scatter¬ 
ing. The mode conversion turns out to occur near specific 
lines in the complex z plane that are called Stokes lines 
in the mathematical literature |28|. A Stokes line is de¬ 
fined as the line where the difference of the WKB phases, 
/ (jfei — k 2 )dC >1 is purely imaginary. Each of the three 
turning points ( p3| ) is origin of three Stokes lines, see the 
Figure. Where the WKB-phase difference is purely imag- 









4 


inary, one of the two modes connected by each turning 
point is exponentially larger than the other. The smaller 
mode cannot be resolved within the WKB approximation 
and may gain a component from the larger mode. The 
single-valuedness of the mode function after a complete 
circle around the turning points uniquely determines the 
conversion rules It follows that the exponen¬ 

tially smaller mode always gains a component from the 
larger one, if the larger mode is present. Therefore, to 
avoid unwanted mode conversion the exponentially larger 
mode should be absent. The Figure shows that the el¬ 
ementary excitations of the sonic white hole can indeed 
remain on the acoustic branch (p5|), which justifies the 
assumptions made to derive the result (jlj). Note that the 
localized acoustic excitations with positive 7 /a represent 
decaying modes for the white hole when a is negative. 
The corresponding growing modes must consist of trans¬ 
acoustic excitations with asymptotics (pif). However, as a 
consequence of the four-fold frequency symmetry of the 
Bogoliubov modes [[To| |24), the spectrum ( 0 ) is deter¬ 
mined by the acoustic modes of the white hole. In con¬ 
trast, the unstable excitations of sonic black holes, if any, 
cannot possess the asymptotics (25) on one of the com¬ 
plex half planes. (24|. Otherwise, sonic black holes would 
be always unstable, and Laval nozzles would be 

unable to stabilize fluids that turn from subsonic to su¬ 
personic speed. Apparently, the opposite process, slow¬ 
ing down supersonic condensates to subsonic speed to 
form a white hole, is intrinsically unstable, generating 
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breakdown shocks 1 25 1 . 

In conclusion, white holes are best avoided in future 
experiments to demonstrate Hawking sound in Bose- 
Einstein condensates. On the other hand, one could still 
use a toroidal geometry [|u|, as long as the elementary 
excitations of the torus do not match the resonances (!li). 
Here the condensate should flow through a constriction 
where it exceeds the speed of sound, establishing a black- 
hole horizon followed by a white hole |uj. Because of 
the periodic boundary condition the spectrum of excita¬ 
tions is restricted. Our theory seems to explain, at least 
qualitatively, why this toroidal arrangement exhibits 
instabilities at well-defined lines in the parameter space 
used. Instead of employing a torus, one could simply 
push a condensate through the optical equivalent of the 
Laval nozzle jl2| , ^5f and let the supersonic quantum gas 
expand into space, like the solar wind |3C|j. Otherwise, 
white holes are as unstable as wormholes [|3l| . 
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